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Abstract 

We study a scalar integro-differential conservation law. The equation was first derived 
in [2] as the slow erosion limit of granular flow. Considering a set of more general erosion 
functions, we study the initial boundary value problem for which one can not adapt the 
standard theory of conservation laws. We construct approximate solutions with a fractional 
step method, by recomputing the integral term at each time step. A-priori L°° bounds 
and BV estimates yield convergence and global existence of BV solutions. Furthermore, we 
present a well-posedness analysis, showing that the solutions are stable in L 1 with respect 
to the initial data. 



1 Introduction 

We consider the initial boundary value problem for the scalar integro-differential equation 

qt+ ( ex p{j^ /(?)) =0, * > 0, *<0, (1.1) 

with initial condition 

q(0,x)=q(x), x<0. (1.2) 

Note that the flux includes a non-local integral term. For notational convenience, we intro- 
duce 



K(q(t,-))(x) = expM°f(q(t,0)d{i 



;i.3) 



The function / : (-l,+oo) — > R € C 2 (R) is called the erosion function. The following 
assumptions apply to /: 

/(0) = 0, f>0, f"<0, lirn /(<?) = -oo, lim ^ = 0. (1.4) 

g->— 1 q^+oo q 
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We remark that the characteristic speed of is 




x = f'(q)K. 

By (|1.3|) and f j 1 . 4 j) . the characteristic speed is always positive, therefore no boundary condition 
is assigned at x = for f)l . 1 f) . 

The equation (jl.ip arises as the s/ow erosion limit in a model of granular flow, studied in [2], 
with a specific erosion function 

Note that this function satisfies all the assumptions in (|1.4|) . In more details, let h be the height 
of the moving layer, and p be the slope of the standing profile. Assuming p > 0, the following 
2x2 system of balance laws was proposed in [12] 

= ^" 1)/l ' (1.6) 
= 0. V ' 

This model describes the following phenomenon. The material is divided in two parts: a moving 
layer with height h on top and a standing layer with slope p > at the bottom. The moving 
layer slides downhill with speed p. If the slope p = 1 (the critical slope), the moving layer passes 
through without interaction with the standing layer. If the slope p > 1, then grains initially at 
rest are hit by rolling grains of the moving layer and start moving as well. Hence the moving 
layer gets bigger. On the other hand, if p < 1, grains which are rolling can be deposited on the 
bed. Hence the moving layer becomes smaller. 

In the slow erosion limit as H/iHl 00 —> 0, we proved in [2] that the solution for the slope p in 
(jl.6p provides the weak solution of the following scalar integro-differential equation 

Pn + exp / — dy\ = . 

V P Jx p(n,y) J x 

Here, the new time variable \i accounts for the total mass of granular material being poured 
downhill. Introducing q = p — 1 and writing t for fi, we obtain the equation (jl.ip with f)1.5|) . 

The result in [2] provides the existence of entropy weak solutions to the initial boundary value 
problem (jl.ip with / given in (jl.5p for finite "time" (which is actually finite total mass). However, 
well-posedness property was left open due to the technical difficulties caused by the non-local 
term in the flux. Furthermore, due to the discontinuities in q, the function k(t, x) = K(q(t, -))(x) 
is only Lipschitz continuous in its variables, therefore one can not apply directly previous results. 
Indeed, classical results as |15j require more smoothness on the coefficients; see also [9]. Some 
closer results can be found in [U [16] where the coefficient k = k(x) does not depend on time. 

In this paper we consider a class of more general erosion functions / that satisfy the assump- 
tions in (jl.4p . and we study existence and well-posedness of BV solutions of (jl.ip . Assuming 
that the slope is always positive, i.e., q > —1, we seek BV solutions with bounded total mass. 
Therefore, we define V = T>c , KQ as the set of functions that satisfy 

V C(h K = lq(x) : inf q(x) > k > -1 , TV {q} < C , ||g|| L i (R _) < C | . (1.7) 



2 



Assume that the initial data satisfies q 6 T>c K0 for some constants Co > 0, kq > — 1. A natural 
definition of entropy weak solution is given below. 

Definition 1 Let T > 0. A function q is an entropy weak solution to (II. li on [0, T] x M_ 

wt/i initial condition ()1 .2j) . if the following holds. 

(HI) g : [0,7] -»• L 1 (M_)n J By(M_), inf^g^x) > -1, and ffce map [0,T] 3 t ^ g(t) is Lipschitz 
in L 1 (M_); 

(H2) q is a weak entropy solution of the scalar conservation law 

f q t + (k(t,x)f(q)) x = 0, 



(li 
(0,x) = q(x) 



with k defined by 



k(t,x) = K(q(t,-))(x) = exp|^°/( 9 (t,0)^} • (1-9) 
Notice that, thanks to (HI), the coefficient k(t,x) in (jl.9p is Lipschitz continuous on [0,T] x 
Now we state the main result of this paper. 



Theorem 1 Assume (|1.4p and let Co > 0, kq > —1 be given constants. Then for any ini- 
tial data q 6 T>c ,k there exists an entropy weak solution q(t,x) to the initial-boundary value 
problem (|l.ip ~ (|1.2p for allt > 0. Moreover, consider two solutions qi(t, •), quit, •) of the integro- 
differential equation M.l}) . corresponding to the initial data 

gi(0,x) = q\{x) , 02(0, x) = g 2 (x) , x<0, 

with q~i, q~2 G 2?c ,ko - ^ en / or an 2/ 7" > t/iere exists L = L(T, Co, «o) > suc/t i/iai 

||«l(t,-)-92(*.-)llLi(R_) < e Lt ||gi -g 2 || L i(R_), ie[0,T]. (1.10) 

Recalling that q = p — 1 = u x — 1, the solution g established by Theorem [1] allows us to 
recover the profile u of the standing layer: 



u(t,x)-x = / q(t,y)dy. (1.11) 

J —00 

Moreover, since JiT x = —Kf(q(t,x)), the equation (jl.lj) can be rewritten as 

Qt ~ K xx = . 

Integrating in space on (—00,2;), using (jl.lip and that K x (q(t, •)) € L^R-.), we arrive at 

u t -K x = u t + Kf(u x -l) = 0. 
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This nonlocal Hamilton- Jacobi equation is studied in [17] , with a different class of erosion func- 
tions /. Assuming more erosion for large slope, i.e., lim g _j. +00 f'(q) = rjo > 0, the slope u x of 
the standing layer would blowup, leading to jumps in the standing profile u. Notice that, in our 
case, only upward jumps in u x can occur as singularities, which corresponds to convex kinks in 
the profile u. 

About the continuous dependence notice that, when k is a prescribed coefficient, the L 1 
stability estimate ([l.lOp holds with L = 0, see (|2.3p . On the other hand, for the integral 
equation (jl.ip . one cannot expect L = in general. Indeed, a small variation in the L 1 norm 
of the initial data may cause a variation in the global term and then in the overall solution. 
However, a special case in which (|1.1U|) holds with L = is when q2 = 0, which indeed is a 
solution of (jl.ip . 

Other problems involving a nonlocal term in the flux have been considered in [lOj [8]. 
Well-known integro-differential equations which lead to blow up of the gradients include the 
Camassa-Holm equation [6] and the variational wave equation [5]. The Cauchy problem for 
(jl.ip with initial data with bounded support is studied in [3] where we use piecewise constant 
approximation generated by front tracing and obtain similar results. 

The rest of the paper is structured in the following way. As a step toward the final result, 
in Section 2 we study the existence and well-posedness of the scalar equation (|1.8p for a given 
coefficient k(t, x). Here k(t, x) is a local term, and preserves the properties of the global integral 
term. Such equation does not fall directly within the classical framework of [15J, where more 
regularity on the coefficients is required (C 1 ). In particular, BV estimates for solutions of (jl.8p 
are needed to obtain the continuous dependence on the initial data, see (I2.23p . We employ a 
fractional step argument to deal with the time dependence of k, and then follow an approach 
similar to [4] (see also where the authors deal with the case of k = k{x) G L°°. We 

further refer to [9] on total variation estimates for general scalar balance laws: their result, in 
our context, would require more regularity (C 1 ) on the coefficient k. 

The properties of the integral operator K, defined at f[ 1 . 3 [) . are summarized in the last 
Appendix. 

2 Local well-posedness of solutions with a given coefficient k 

In this section we study the well-posedness of the scalar equation (11.8P for a given coefficient 
k(t,x), by reviewing some related results and completing the arguments where needed. 
Throughout this section, we will use u as the unknown variable. Consider 




t > 



(2.1) 
(2.2) 



u(0, x) = u(x) 



x < 



where k(t,x) satisfies the following assumptions, for some T > 0: 

k(t,x) € L°° ([0,T] x R_), it is Lipschitz continuous and inf 4jX k > 0; 
(K) TV {k(t, •)}, TV {k x (t, •)} are bounded uniformly in time; 
[0,T] 9 t -> k x (t, ■) S L 1 (M_) is Lipschitz continuous. 
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The above assumptions on k are motivated by the properties of the integral operator K, see 
Proposition [2] in the Appendix. 



Theorem 2 Assume f satisfies (II. 4h and k(t,x) satisfies (K). Let Cq > 0, kq > — 1 be given 
constants. Then there exist two constants C\ and K\, with possibly C\ > Co and —1<hi<kq, 
and an operator P : [0, T] x T>c ,k ~ * ^Ci,« a such that: 

1) the function u(t,x) = Pt{u) is a weak entropy solution of (|2.ip with initial data u(0, •) = 
u e T>c , K0 ; 

2) for any u\, Ui G T>c 0)K0 one has 

\\Pt(ui) - Pt(u2)\\hiQi_) < \\ui -n 2 ||Li(M_) • (2-3) 

Proof. Let u S T>c ,n - We introduce the parameter At > and define t n = nAt for any 
integer n > 0. We approximate the coefficient k by 

kAt(t,x) = k(t n ,x) (t,x) G [t n ,t n +x) x R_ , n>0 (2.4) 

which is constant in time on each interval [i n ,i n +i), and consider the equation 

u t + (k A t(t,x)f(u?) =0. (2.5) 

By adapting the analysis in [3], on each interval [i n ,i n +i) the entropy solution for (j2.5j) . call it 
UAt, exists and the corresponding operator (t,u) i-> UAt(t, •) is contractive in L 1 ^.), provided 
that UAt is bounded from both below and above. Furthermore, the complete flux 

F(t,x) = kAt(t,x)f(uAt(t,x)) 

has the following properties: its sup norm does not increase in time, 

\F{t,x)\ < sup|F(t„,-)|, te{t n ,t n+1 ), (2.6) 

as well as its total variation: 

TV{F(t,-)} < TV{F{t n ,-)}, te{t n ,t n+1 ). (2.7) 

We now establish the lower and upper bounds for UAt- F° r notation simplicity, in the 
following we denote by k(t, x) and u(t, x) the approximate coefficient and solution respectively, 
without causing confusion. We define the constants ko, L, L\ such that, recalling (K), one has: 

/c = infA;>0; (2.8) 

t,x 

\k(ti,xi) - k(t2,x 2 )\ < L (\ti - t 2 \ + \xi - x 2 \) for all ij, x i} i = l, 2; (2.9) 
TV{fc(ti, •) - k(t 2 , ■)} = \\k x (h, •) - k x (t 2 , .)|| L i(R_) < L i\h ~ h\ (2.10) 



5 



and set L2 = L/ko- We first give some formal arguments. The evolution of the complete flux 
F = kf(u) along the characteristic x(t) with x = f'(u)k follows the equation 

j t F(t, x(t)) = (kf) t + f'k(kf) x = ktf = jF . (2.11) 

By our assumptions (K), the term kt/k is uniformly bounded. Therefore, \F\ grows at most 
at an exponential rate, and remains bounded for finite time t < T. Therefore |/(it)| remains 
bounded as well. By the 4 th assumption in (jl.4p . it never reaches —1 in finite time, leading to a 
lower bound on u. 

The same argument leads to an upper bound for /(it), if /(it) — > +00 as u — >■ +00. However, 
if /(it) — > fo > as u — > +00, we need a different argument. We observe that, along a 
characteristic x(t), one has 

^-u(t,x(t)) = -k x (t,x)f(u). (2.12) 
at 

By the lower bound on it, the growth of it remains uniformly bounded, yielding an upper bound. 
We now make these arguments rigorous for the approximate solutions. At time t = one 

has 

\k(0,x)f(u(x))\ < Ci (2.13) 

for some C\ > that depends on the bounds for k and u. We claim that, as long as the 
approximate solution exists, we have 

\k(t,x)f(u(t,x))\ < de L2t . (2.14) 

Indeed, by (I2.13D and (|2.6h . the inequality (12.141) is valid on [0, t\). Assume now that (|2.14j) is 
valid on [0, t n+ \), n > 0, i.e., 

\F(t,x)\ = \k(t n ,x)f(u(t n ,x))\ < Cie L2 *", te[t n ,t n+l ). (2.15) 

At time t = one has 

\k(t n+1 ,x)f(u(t n+1 ,x))\ = k ^". +1,X ? \k(t n ,x)f(u(t n+1 ,x))\ 

K[t n , X) 

< fl + -^AtVsup|fc(t n ,rE)/(n(t n+1 ,x))| < e i2A * • C x = C ie L ^ 

\ k J x 

By induction, this proves (I2.14p . which in turn gives the lower bound k\ for it. The upper bound 
also follows if /(it) — > +00 as it — > +00. 

Finally, we consider the case that /(u) — > fo > as it — > +00. At any given point (t, x) one 
can trace back along an extremal backward generalized characteristic, which is classical on each 
(t n ,t n+ i) and continuous up to t = 0. Since now the r.h.s. of (|2.12p is bounded, then u grows 
at a linear rate, and therefore remains bounded. 

We remark that the lower bound yields an a-priori bound on the wave speed. Indeed, since 
/' is a decreasing function, the characteristic speed is bounded, 

A = *;/'(«) < HfcHoo/'fci). 
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Bound on total variation. We estimate the total variation of F(t,x) = k(t,x)f(u(t,x)). On 
the interval (t n ,t n+ i) the coefficient k is constant in time and we use ()2.7p . On the other hand, 
the total variation might increase at t n when k is updated. Then we observe that 



F(t n ,x) 



1 + 



therefore 



TV{F(t nr )} < 1 + 



k(t n ,x) k(t n —i,x) 

k(t n —l, x) 
\k(t n , •) — k(t n —i, -Jlli 



inf k(t n -i, ■ 



F(t n , x) , 



TY{F(t n -,-)} 



(2.16) 



Thanks to ([2T8jl - ([2TT0jl . we have 

\\k(t n , •) - fcfa-i, Qlloo 
inf k(t n -i, •) 



+ sup \F\ ■ TV 



< L 2 At, TV 



k{t n , •) k(t n —\, •) 

k(t n , •) — k(t n —i, •) 
k(t n -i, •) 



(2.17) 



for a suitable constant L3 independent on At. Moreover F = kf is uniformly bounded thanks 
to (K) and the bounds on u. Hence we conclude that 

TV{F(t n ,-)} < (1 + L 2 At) TV {F(t n _i, •)} + L 4 At 

for a suitable L4 > 0. By induction it follows that 

TV{F(i,-)} < e L2t TV{F(0+,-)} + ^(e L2 '-l) . 

Recalling that f(u) = F/k, one obtains the BV bound for f(u(t)), 

(inf /') TV {«(*,-)} < TV •))} < 4l TV ^-)> + ?Rfe TV -[^-)}- 

mi k (inf 



This gives a bound on the total variation for u(t): 

TV{«(t)} < C[TY{F(t,-)} + TY{k(t,-)}\ < C\(t) 



(2.18) 



where the constant C depends on inf^ u, sup x u, vai x k, sup x k. Hence the total variation of u 
may increase in time but it remains bounded as long as u remains bounded. 

Taking the limit At — > 0, the coefficient kAt converges uniformly to k. Correspondingly, 
the family uj\t converges (up to a subsequence) to a weak solution u of the original equation, 
satisfying the same upper and lower bounds and fj2. 18f) . 

Moreover, in the limit as At — > 0, the Kruzkov entropy inequalities for equation (|2.ip 



d t \u-a\ + d x [k(x,t)\f(u) - f(a)\] + sign(u-a)k x (x,t)f(a) < (2.19) 
for all a£l, hold in the sense of distributions. □ 
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Next we establish the continuous dependence on the coefficient function. We rely on a result 
in |14j (Corollary 3.2) that applies to Cauchy problems and to the case of k = k(x), that is, the 
coefficient does not depend on time. 

For convenience of the reader we report that statement of [2] adapted to our situation. 
Consider the two equations 



u t +[kf(u)) = 0, t>0, (2.20) 

u t +(kf(u)) =0, t>0. (2.21) 



Proposition 1 For i£l, let k(x), k(x) G BV(M) satisfy 

k x , k x G BV(R) ; inf k , inf k > a > 

for some positive a. Consider the initial data uq, uq G BV(M.) for the two equations (j2.20|) . 
(|2.2ip respectively and let u(t,x), u(t,x) be the corresponding solutions, assuming that they are 
bounded from above and bounded away from — 1. Let C\ be a bound on |/| over the range of the 
solutions. Then 

\\u(t, -)-u(t, -)IIli(r) < ||«o - «o||li(R) 

+ t{c 1 TV{k-k} + C 2 (l + TVu + TVu ) ||ib-jb||oo} (2-22) 

where C2 depends on the bounds on u, k, TV {A;} and on u, k, TV {k}. 

The continuous dependence property for our problem follows from Proposition [U by properly 
extending the IBVP into Cauchy problems. 

Theorem 3 For x < 0, let k(t,x), k(t,x) satisfy the assumption (K), and assume that the 
initial data u belongs to T>c ,n (defined at (|1.7p ). Let u(t,-), u(t,-) be the solutions of the 
conservation laws (|2.20p . (|2.2ip respectively, with the same initial data u, for some time interval 
[0,T] (T>0). 

Then, the following estimate holds 



7lKV)-£(V)|| L i(R_) < C l sup TV \k(t, •) — k(t, •) j 
1 te[o,T] 1 ' 

+ C 2 ^l + supTVu(r, •) +supTV«(r, -)^) \\k - fc||L°°([o,t]xR_) , ( 2 -23) 



where C\ is a bound on \ f\ over the range of the solutions and C2 depends on the bounds on the 
solutions, the coefficients and their total variation TV {k(t, ■}, TY {k(t, •}. 

Proof. The IBVP (|2.ip - (|2.2p can be extended to the following Cauchy problem 

ut + (k(t, x)/(«)) = , xeR, t>0, (2.24) 
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with extended initial data 



u(0, x) 



and the extended coefficient function k(t, x) 



u(x) 
u(0- 



k(t, x) 



lim k(t, y) 



for x < . 
for x > 



for x > . 



(2.25) 



Due to the fact that the characteristic speed is positive, the solution for the Cauchy problem 
(pTMl) - (jiI25]) restricted oni<0 will match the solution for the IB VP (f2TT|) . 
In a same way, the IBVP (|2.2ip is extended to the Cauchy problem for 



ut + (k(t,x)f(u) 



0. 



x € 



t > 



(2.26) 



with data (|2.25p . Without causing confusion, let's still denote u(t, x) and u(t, x) the solutions 
for (|2.24p and ()2.26j) . respectively, and let u&(t,x) and u&(t, x) be the corresponding approx- 
imate solutions, constructed in the same way as in the proof of Theorem with approximate 
coefficients kAt and k&t as in (|2.4p . 

Denote the distance between these two solutions by 

eA(t) = ||«a(*,0 -«a(*,-)|| l i(K) • 

Notice that eA(0) = and that e&{t) > ||ua(£, •) — u&(t, •) • 

On each time interval [t n ,t n+ i) the coefficient is constant in time and the assumptions of 
Proposition [TJ are satisfied. Hence, from ()2.22p . we have the following estimate 

eA(Wi) - eA(*n) < At C*iTV R |fcA t (t n , •) - ^At(tn, •)} 

+ At C 2 (l + TV R _n(t n ,-)+TV K _n(t n ,-)) k At (t n , ■) - k At {t n , ■) (2.27) 

L°°(R) 

for some constants C\ and C2 that are uniform on [0, T]. Notice that, in the above lines, 
TVir I A; At — ^At| coincides with TVk_ of the same quantity and, similarly, the L°°-norm on 
R coincides with the L°°-norm on R_. Concerning TVru (similarly for TViS), we replaced it 
with TVi_it with an error that is bounded and possibly depending on T. 
Summing up (|2.27p in n, we get 



N-l 



eA(*iv) - eA(0) 



V" e A (Wi) - e A (t n ) < t N Ci sup TV R _ \ k A t - k A t\ 

te[o,t N ] L J 



n=0 



+ t N C 2 ( 1 + sup TV R_-u(t, •) + supTVM_«(t, 

t t 



kAt ~ kAi 



L°°([0,tjv]xR_) 



Now taking the limit At — > 0, we get (12, 23ft . completing the proof. 



□ 
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3 Well-posedness of the integro-differential equation 



In this section we prove the main Theorem[TJ In Subsection l3.1l we define a family of approximate 
solutions to (|l.ip — (|1.2p and show their compactness, locally in time. Then in Subsection 13.21 we 
show that the limit solution can be prolonged beyond the existence time, by improving the 
estimates on upper and lower bound for the exact solution of (|l.ip - (|1.2p . Finally, in Subsec- 
tion 13.31 we show that the flow generated by the integro-differential equation (jl.ip is Lipschitz 
continuous, restricted to any domain T> given at (II. 7p . 

3.1 Local in time existence of BV solutions 

In this Subsection we prove the following existence theorem. 

Theorem 4 Let Co, kq be given constants and let q(x) G L 1 (1R_) n BV(§&-) such that 
(a) mi x<0 q(x) > k > -I; 



admits an entropy weak solution q(t, x) on [0, T] x R_ that satisfies 

(a) ' ra£ x<0 q(t,x) > «i>-l; 

(b) > TV {<#,-)} <Ci ; 

(c) ' lk(t,')llLl(R_) < ll?llLl(R_) • 

Proof. We define a sequence of approximate solution to the scalar equation (|l.ip - (|1.3p . We 
fix At > and set t n = nAt, n € N. The approximation is generated recursively, as n starts 
from and increases by 1 after each step. For each step with n > 0, let q(t,x) be defined on 
[0, t n ) x R_ and set 



(b) TV {$(•)}< C ; 



( c ) II?IIli(r_) < Co • 



Then there exist T > ; K\ > — 1 and C\ > such that 




(3.1) 




Then we define q on [t n ,t n +i) x M_ as the solution of the problem 
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This procedure leads to a solution operator t i— > S At q = q At (t, •), defined up to a certain time 
T = T(At, q) > 0, of the problem 

f q t +{k A %x)f(q)) x = 0, i>0 
\ q(0,x) = q(x) , 

where k = k At is defined by 

k At (t,x) = ^X[Wn + i)(*)" fc n(x). (3.3) 

n>0 

Notice that the operator S At has the semigroup property S At + = S At S At for t\, t% G (At)N. 
Now we prove uniform bounds, independent of At, on the family of approximate solutions. 

The L l bound. This follows by the application of (|2.3[) in Theorem [2j at each time step 
[t n , t n+ i), and the fact that 1 1— > q(t, •) is continuous in L 1 . Until the solution is defined, we have 

h(t,-)hi < lk(o,-)l| L i- (3.4) 

Lower and upper bound on q. Define 

z(t) = ini q(t,x), w(t) = sup q(t,x). 

X X 

We observe that, by comparison with the equilibrium solution u = 0, (i) if z(0) > then z{t) > 0; 
and (ii) if w(0) < then w(t) < for all t > 0. 

Now consider —1 < z(0) < and w(t) > 0. Choose 5 and M such that z(0) > —1 + 25 and 
w(Q) < M/2. For example, one can take 5 = (k + l)/2 and M = 2w(0). Let T = T(S, M) > 
be the first time that one of the following bounds fails, 

z(t)>-l + S, w(t)<M. (3.5) 

Then, for t < T, from the analysis of equation fl3.2f) (see (]2.12p ). we find that z and w are 
continuous and satisfy 

z(t) > z(0)+sup\k At (t,x)\ [ f(z(r))dr, z<0, (3.6) 



x J o 



w(t) < w(0) +sup\k At (t,x)\ / f(w(T))dT, w>0. (3.7) 





Note that in (|3.6p we have f(z) < 0, and in f)3.T|) we have f(w) > 0. For |fc^*|, we have the 
estimate 

\k A %x)\ = \k At (x)f(q(t n ,x))\ < expU°\f(q(t n ,0)\d^\f(M) 



< f(M)exp{f'(-l + 6)\\q\\ hl } < C(S,M) 
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This gives us 

z(t) > z(0) + C(S,M) [ f{z{j))dr > z(0) - C{5, M)t |/'(-l + S)\ , 

Jo 

w(t) < w(0) + C(5, M) [ f(w(T))dT < w(0) + C(5,M)tf(M). 

Jo 

We conclude that the bounds in (13. 5p hold for t <T with 

T(6,M) = min{T 1 ,r 2 }, 

where 

Tl(5 ' M) = C(5,M)\r(-l + 5)\ ' T2(5 ' M) = C(6,M) 2 f(M) ' 
yielding the lower and upper bounds. 

Finally, if z(0) > and w(0) > 0, or if z(0) < and w(Q) < 0, then we would only need to 
establish one of the bounds in (|3.5j) . and the result follows. 

Bounds on f,f',k. Once we have a lower, upper bound on q and the bound on H^Hl 1 ) we 
immediately find that 

f(q(t,x)), f'(q(t,x)), I f(q(t,0)<% G L°° ([0, T] x R_) (3.8) 



uniformly w.r.t. At. By definition f|3.3|) of /c, we can easily verify that the following properties 
hold uniformly w.r.t. At: 

(i) k £ L°° ([0,T] x R_), inf t>af fe > 0; 

(ii) ft, €L°°([0,r] xl_); 

(iii) TVfc(t, •) is bounded uniformly in time . 

Indeed, (i) follows from the definition of k and (|3.8[) . About (ii), at each time t we have 
k(t, ■) = k n {-) for some n, and for = —k n f(q(t n ,-)). Then /c x € L°° because of (i) and (|3.8p . 
Finally 

TVfc(t,-) = H^IIli = ||fcn/(g(*n,-))llLi < M \\k\\oo \\q(t n , < M HfclloollgllLi 

where M = sup /', that depends on the lower bound on q. 

Lastly, from (i) and (|3.8p one obtains a uniform bound on the characteristic speed kf'(q). 

Bound on the total variation of q. By definition of the total variation 

f-0 



TV{q(t,-)} = lim \ ! \q(t,x)-q(t,x-h)\ 
h^0+ h J_ 00 



dx . 
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we have, for any h > 



h 



l(t,x)-q(t,x-h)\dx < TV{q(t, ■)} . 



(3.9) 



The total variation of q does not change at time t n when k is updated. Now consider a time 
interval t E [t n , tn+l); an d we estimate the change of the total variation of q in this time interval. 
We have 

rO rO 

\q(t n +i,x) - q(t n +i,x - h)\ dx < / \q{t n , x) - q(t n , x - h)\ dx 

) J — oo 

/•tn+l 

+ / £(r)dT 



(3.10) 



where 



S(t) = limsup 
6» — s-0-1- 



J-oo \q( T + 0,x - h) - q(r + 9, x)\ dx 



Here q is the entropy solution to 

( u t + (k n (x)f(u)) x = 0, t>T,X<0 
I u(t, x) = q(r, x — h) . 

On the other hand, q(t, x — h) is a solution of 

\ u t + (k n (x - h)f{u)) x = , t > t , x < 
| u(t, x) = g(r, x — h) . 

Using the estimate (I2.22P we find 

S(t) < \\f\\ O0 TV{k n (--h)-k n (.)} + C(l + TV{q(T,-)})\\k n (--h)-k n (- 

for a suitable constant C. Notice that 



(k n )x(r, y) dy 



x—h 



< h\\k n f\\ c 



\k n (x - h) - k n (x)\ = 
and that 

TY {k n (- - h) - k n (-)} = \\(k n ) x (--h)-(k n ) x (-)\y 

< mn{ . _ h) _ fc n (.)) / ( g ( T) .))|| L1 + ||^(- -h)-(f (q(r, --h))-f (q(r, .)))|| L a 

< /ill^n/Hoo • ||/(9(T, -))IIli + ll fc n|| L °o ll/'lloo ||?(T, •) -^(T, ' -/l)|| L i • 

In conclusion, using also H3.9j) . we obtain 

£(r) < ^|M 1 + M 2 TV{ g (r,-)} + M 3 i||g(T,-)-g(T,--/ i )|| Ll | 
< h {Mi + (M 2 + Ms) TV {g(r , •)}} 
where Mj depend only on a-priori bounded quantities. Now from (I3.10P we obtain 



TV{q(t n+1 ,-)} < TV{g(i n ,-)} + 



Mi + (M 2 + M 3 ) TV {q(r, •)} dr. (3.11) 



We conclude that the total variation of q may grow exponentially in t on each interval (t n , t n +i) 
but it remains bounded for any bounded time t. 
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Convergence to weak solutions; Existence of BV solutions. Now, without causing 
confusion, we will use q A (t,x) for the approximate solution, where A = At is the step size. Let 
k A be the approximated coefficient of the equation, defined in (|3.3[) . 

By compactness, there exists a subsequence of {q A (t,x)}, as A — > 0, that converges to a 
limit function q(t,x) in L[ oc . Let k(t,x) be the integral term, (jl.9p . corresponding to q, which 
is uniformly bounded as well as the k A . We have 

k A (t, x) - k(t, x) = 0(1) { f f(q A (t n , £)) di - f f(q(t, 0) di 

v. J X J X 

= 0(1) jsup TV {f(q A (r, •))} sup ±A + f\f{<f& 0) ~ /(<?(*, 0)] # 

that vanishes as A — > 0. Therefore we can pass to the limit in the weak formulation. On the 
interval [i n ,i n+ i], q A (t,x) satisfies 

/ / (q A 4>t + k A f(q A )(p x )dxdt= [q A (p(t n+1 ,x) - q A cp(t n ,x)} dx 

J t n J — oo J — oo 

for some test function <p with compact support inside [0,T] x M_. Summing this up over n, we 
get 

T r r0 

{q A (t ) t + k A f{q A )^ x )dxdt = / [q A (f)(T, x) — q A (j)(0, x)] dx . (3.12) 



o 

Since q A — >■ g in L[ , /((? A ) — >■ f(q) in L^ oc , /c A — ^ pointwise and fe A , k are uniformly bounded, 
by dominated convergence we can take the limit as A — > and have the convergence of (|3.12p 
to 

J J [q(t,x)4> t (t,x) + k(t,x)f(q(t,x))4> x (t,x)]dxdt = / [qcj)(T,x) - q<f)(Q, x)} dx . 

JO J —oo J —oo 

This completes the proof of existence of BV solutions for (jl.ip . □ 
3.2 Global existence of BV solutions 

Once the BV solutions exist locally in time, we can further show that they enjoy better properties 
than the ones deduced from the approximate solutions. In particular we show that the lower and 
upper bounds on q do not depend on time t, leading to global in time existence of BV solutions. 

Let q be an entropy weak solution of (|3.ip on [0, T] x M_. We will now improve the needed 
bounds. 

Lower bound on q. Given any point (i,x) G (0, T) x R_, let t — > x(t) be the minimal 
backward characteristic (which is classical), defined for t G [0,t\. By setting q(t) = q(t,x(t)), we 
have 

f x'(t) = k(t,x)f'(q{t)), x(t)=X, 

{ (3.13) 
{ q'(t) = -k x (t,x(t))f(q) = kf(q) 2 >0, q(t) = q(t,x-) . 

We see that the solution q is non-decreasing along any characteristics. Therefore, we have 
inf;; q(t, x) > inf q(x) > kq > — 1 for all t > 0. 
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Upper bound on q. Again, consider a point (t, x) and let t — > x(t) be the minimal backward 
characteristic through it. From the second equation in f)3. 13j) we see that if g(0,t(0)) < 0, then 
q — > as t — > +00. Now consider q(0,x(0)) > 0, and we have q(t,x(t)) > for all t. Define 

W(t,x)= ! \q(t,y)\dy, x<0, (3.14) 
J —00 

that satisfies 

< W(t,x) < ||?(0,.)||li(r_). 
Using (|2.19p with a = 1, we have 

/X px 
\q(t,y)\ t dy < - / (k(t,x)\f(q)\) dy = -k\f(q)\. 
-00 J -co x 

The variation of W along the characteristic is 

jW(t,x(t)) = W t + x'W x < -k\f\ + \q\kf = k{-\f\ + \q\f'(q)) 

. i( _ /+ ^ )) ._^pa._-(^ fJ ). ( , 15) 

Here we remove the absolute value signs because q > 0. Then, (|3. 15[) implies that 

WMt)) + lWMfm = c 



along characteristics. This gives the bound 
q(t,x(t)) q(0,x(0)) 



+ W{0,x(0))-W(t,x(t))<C 1 , (3.16) 



f(q(t,x(t))) /( 9 (0,x(0))) 
where C\ can be chosen independently of (i, x). Recalling (|1.4p . we have 

lim = +00 . 

Therefore, fj3. 161) implies an upper bound for q for all t. The uniform bound on the total variation 
follows because the constants Mi in (|3.11|) are now bounded uniformly in time. 

3.3 Continuous dependence from the data for the integro-differential equa- 
tion 

In this section we prove the last part of Theorem [TJ showing that the flow generated by the 
integro-differential equation (|1.1|) is Lipschitz continuous, restricted to any domain T> C L 1 (M_) 
of functions q(-) satisfying the following uniform bounds in (jl.7p . for some constants Co, n$. 

Consider two solutions q±(t, •), q%(t, •) of the integro-differential equation (jl.ip . say with 
initial data 

qi{0,x) = qi(x) , q 2 (0,x) = q 2 (x) x<0, 
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and satisfying the conditions in (|1.7p for t E [0, T]. We are going to prove that 



ki(t,0-«2(v)ll 



+ L 



\qi(s, •) - q2(s, 



ds . 



(3.17) 



for a suitable constant L. By Gronwall lemma, this yields (jl.lOp . hence the Lipschitz continuous 
dependence of solutions of (jl.ip on the initial data. 

Define the functions k\(t,x), k2(t, x) as in (jl.9p . corresponding to qi(t, x), q2(t,x) respec- 
tively. Now set 

f h(t,x) if te [0,0], 

{ k 2 (t,x) if t>9. 

Finally, for any given 6 E [0, T], let q 9 = q e (t,x) be the solution of the conservation law 

q t +(k e (t,x)f(qj) =0, q 9 (0,x)=q 2 (x). (3.18) 



Observe that, for each fixed 9, the distance between any two entropy-admissible solutions of the 
conservation law (|3.18p is non-increasing in time. In particular, for 9 = T, call q the solution of 

g t + (fei(t J a;)/(g)) !e = 0, 

with initial data q(0,x) = q2{%) (see Figured]). We have 



\vi(t, -)-q(t, Oil 



< 



m ~ <?2|| L i (R _) 



for all t E [0, T] . 



(3.19) 











! 2 








q° 


^1 


















— " ( l 






; Mj 



Figure 1: The flow of solutions qi,q,q 9 ,q2 for the integro-differential equation. 

Moreover we can use the Lipschitz property of the solution operator for (jl.lj) with k = k 2 
fixed, and get the distance estimate 



\q(T,-) -«j(r,OllLi(tt_) 
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< 



£7(t) dr 



(3.20) 



where 

E(r) = limsup 11^ + M - $(r + Mil., 

Indeed, observe that q(r, •) = q e (r, •) whenever r < 9, for any r € [0,T]. 

To compute the integrand in (13.20p . observe that the functions h i— > q T (r + /i, •) and h i— > 
g(r + /i, •) take the same value q(r, •) when /i = 0, and h ^ q T {t + h, x) satisfies the conservation 
law 

q h +(k 2 (T + h,x)f(q)) x =0, (3.21) 

while h i — y q(r + h, x) solves 

q h +(k 1 (T + h,x)f(q)) x =0, (3.22) 

for h > 0. By using ()2.23|) in Theorem [3j we can measure the error term E(t). By the facts 
that ||<f (r, •) || LOO , ||g(r, •) H^, TV {f(r, •)}, TV {q(r, •)}, TV {h(r, •)} and TV {fc 2 (r, •)} are all 
bounded, the coefficients C\ and C 2 in ()2.23j) are all bounded constants. Let M be a generic 
bounded constant, we get 



sup TV (k^t, •) - fc 2 (V)) + \\h - Ml°°([t,t+Mx]R_) 

T<t<T+h 



\h(T, ■) - k 2 (r, ^Hloc = M-sup 

x<0 



f0 

- / / (Q2(r,y)) dy 



\\q T (T + h r ) -g( r + /i,-)|| L i < M/i 
Therefore, we have 

£(r) = M • TV {^(r, •) - A; 2 (t, ■)} + M ■ ^(r, •) - k 2 (r, .)|| Loo (3.23) 
Recalling the definitions of fej, £; 2 we deduce that 

f°fMr,y)) 

J X 

= M-\\q 1 (T,-)-q 2 (T,-)\\- L i, (3.24) 

and, using also (13,24p . 

TV{fei(r,.)-fc2(r,-)} = ||(fcl)*(r,-)-(AHj)x(r,-)||Li 
= llfci(^r, ■) / («i(t, *)) - *a(-r, -) / (^(t-, 

< \\(h(r, •) - k 2 (r, •)) / (gi(r, -))|| L i + |Mr, •) • (/ (gi(r, •)) - / Mr, 0))|| L i 

= ll^i (t,0 - fair, oHloo • ||/(9i(t,-))|| L i + ll fc 2(T, -)II L - lki( T >-) -92(r, -)ll L i 

= M|| gi (r,-)-g 2 (r,-)|| L i • (3.25) 
Putting the estimates (|3.24p and (|3.25|) into (|3.23j) . we get 

E(r) < L-\\ qi (T,-)-q 2 (T,-)\\v 
for a suitable constant L. Inserting this estimate in (I3.20P and using (I3.19P one finally obtains 
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A Properties of the integral operator 



In this Appendix we prove some properties of the integral term k in terms of a Lipschitz flow 
1 1 y q(t, •). The operator K, see (|1.3p . is defined on the set 



DBV(l 



inf q(x) > — 1 
x<0 



and valued in Lzp(R_). Its properties are summarized in the following Proposition. 

Proposition 2 Let Co, kq, T be given positive constants. Assume that the map q : [0, T] 
Dc ,k is Lipschitz continuous as a function in L 1 (M_). 
Define k as in (|1.9p . Then 



(K) 



k(t, x) : [0, T] x M_ — > R + is bounded and Lipschitz continuous, with inf^ A; > ; 
TV k(t, •) , TV k x (t, •) are bounded uniformly in time; 
[0,T] 9 i -> fec(t, •) € L 1 (M_) is Lipschitz continuous. 



Proof. To begin, notice that the quantity k is well-defined and is Lipschitz continuous on 
[0,T] x M_. 

Let L be a Lipschitz constant of the map [0,T] 9 t q(t) G L 1 (R_). From the bounds ([LTD 
one easily deduces that 



||9(*j*)||l°°(r_) < Co, 

||/(?(<,0)lk»(R_) < max{|/(C )|,|/(« )|} , 

II/(?(*,-))IIli(r_) < I/WI-II^OIIlhk^^CoI/'MI, 

ll/(?(*l,0) -/(9(*2,0)IIli(M_) < £|/'(«0)|-I*l-*3|- 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



By the assumptions on q we find that 

r0 



< \\f(q{t,-))\\ 



< C |/'(«o)|. 



Hence the integral term k is bounded and satisfies 

< exp (-C \f'(n )\) < k(t,x) < exp (C7 |/'(«o)|) 



Moreover, for all < t\ < t 2 we have 
r 



[f(q{h,0)-fi<l(h,0)) 



< Wf(q(h,-))-f(q(t2,-))hH^) < L\f'(K )\-\h-t 2 \. 



This leads to the Lipschitz continuity in t for k(t,x). Namely, for all x we have 



\k(t u x)-k(t 2 ,x)\ = 0(1) 



[/(9(*i,o)-/(?(*2,o)] <*e 



<L |ti-t 2 |. 



(A.5) 



Here the Lipschitz constant L depends on the parameters L, Co, «o- 
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From the definition of k, the derivative function k x satisfies 

k x = -kf(q) G L 1 n L°° . (A.6) 

This immediately shows three facts: (i) k(t,x) is Lipschitz in space variable x, (ii) k(t,-) € 
BV(M_) where the BV bounds are uniform in t, and (iii) k x (t, •) € BV(M.—). 
From (|A.6j) we get the estimate on the total variation of k x 

TV (k x ) < TV (k) • II/^Hloo^.) + |N|l~(imTV (f(q(t, •))) < M TV (q) , 

with M depending on the parameters L,Cq,ko- 

Finally, we show that [0, T] k x (t, •) 6 L 1 (1R_) is Lipschitz continuous. By using (|A.5|) . 

fOj) and (fOj) . one has 

H*«(ti,-)-Ax(t2,0llLi(R_) = TV{fc(ti,-)-fc(t2,0} 

= •)) " fc(«2, 0/(?(t2, -))IIl1(R_) 

< •) - fc(i 2 , -)lk<-(K_)II/(^(*i, 0)IIlhr-) 

+ WHh, -)IIl<»cr-)II/(?(<i. •)) - /(<z(*2 S -))h H ^) 

< M\t 1 -t 2 \ 

with M depending on the parameters L, Cq, kq. □ 
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